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We define the non-perturbative part of a quantity as the difference between its numerical value 
and the perturbative series truncated by dropping the order of minimal contribution and the higher 
orders. For the anharmonic oscillator, the double-well potential and the single plaquette gauge 
theory, the non-perturbative part can be parametrized as A e~^^^ and the coefficients can be 
calculated analytically. For lattice QCD in the quenched approximation, the perturbative series 
for the average plaquette is dominated at low order by a singularity in the complex coupling plane 
and the asymptotic behavior can only be reached by using extrapolations of the existing series. 
We discuss two extrapolations that provide a consistent description of the series up to order 20-25. 
These extrapolations favor the idea that the non-perturbative part scales like (a/ro)^ with a/ro 
defined with the force method. We discuss the large uncertainties associated with this statement. 
We propose a parametrization of ln((a/ro)) as the two-loop universal terms plus a constant and 
exponential corrections. These corrections are consistent with a1_i„^p and play an important role 
when /3 < 6. We briefly discuss the possibility of calculating them semi-classically at large (3. 

PACS numbers: 11.15.-q, 11. 15. Ha, 11.15.Me, 12.38.Cy 



I. INTRODUCTION 

The perturbative renormalizability of the standard 
model is a remarkable property that played an important 
role in establishing its phenomenological prominence. It 
guarantees that we can calculate the radiative corrections 
at any order in perturbation theory. However, Dyson's 
argument 0] suggests that these perturbative series are 
divergent and one needs to truncate them in order to 
get a finite result. A procedure often used for this pur- 
pose consists in dropping the smallest contribution at a 
given coupling and all the higher order terms. The dif- 
ference between this truncated series and the numerical 
value of the quantity calculated can be called the non- 
perturbative part (NPP). The NNP depends on the trun- 
cation procedure and also on the renormalization scheme. 
However, if the NPP has a simple parametrization and if 
we can calculate it approximately, we will have made a 
big step toward a complete solution of the problem. 

In this article, we discuss the NPP as defined above 
for a variety of models. Our main goal is to find an ac- 
curate formula for the average plaquette in lattice QCD 
in the quenched approximation. The paper is organized 
as follows. In Sec. ^ we calculate the minimal error 
that can be made at a given coupling, if we assume that 
for sufficiently small coupling, the error is dominated by 
the first order dropped. In Sec. IIIII we show that the 
general estimates of Sec. m work well for the anhar- 
monic oscillator. In Sec. IIVI we show that for the ground 
state of the double-well potential, the one instanton ef- 
fect is larger than the error estimated on the basis of the 
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asymptotic behavior of the perturbative series. However, 
for the average of the two lowest energy levels, we ob- 
tain a good agreement with the general estimate. The 
case of a SU (2) lattice gauge theory with one plaquette 
is discussed in Sec. In that case, the NPP can be 
calculated exactly and it can be checked that the error 
formula works well. In contrast to what happens for the 
anharmonic oscillator where the error at a given order 
has a given sign, the examples discussed in Sees. IIVI and 
13 show more complicated sign patterns that can studied 
easily because we can calculate the numerical values very 
accurately. In particular. Figs. ^ to IHl should be under- 
stood as a tutorial for the interpretation of Figs. 1141 and 

In all the above examples the NPP can be approxi- 
mated as A e~'-"/^ with A a generic notation for the 
expansion parameter. Could the same type of result be 
obtained in lattice gauge theory? Lattice QCD in the 
quenched approximation is introduced in Sec. IVll Sim- 
ple hypotheses (power of the two loop renormalization 
invariant scale, power of the lattice spacing) for the NPP 
of the plaquette are compared with the numerical data 
using a series calculated up to order 10 This order 
is not large enough to decide if any of the hypothesis 
is adequate. Higher order extrapolations are necessary 
to discriminate among the various possibilities. In Sec. 
IVIII two extrapolations are discussed. One is based on 
the existence of complex singularities jaL the other on 
the existence of an infrared renormalon 0, |^ 0] . These 
two extrapolations predict quite correctly the 16-th co- 
efficient as it can be obtained from a graph in Ref. 0- 
They both seem consistent with the possibility 0,13 that 
the NPP is proportional to the 4-th power of the lattice 
spacing as parametrized using the force method in Refs. 
[3, ^3 ■ The uncertainties associated with the extrapola- 
tions are discussed in Sec. IVIIII In Sec. IIXI we discuss 
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an exponential parametrization of the lattice spacing as a 
function of the inverse coupling. Our result is consistent 
with the possibility of corrections scaling approximately 
like a^_/™,, in good agreement with a suggestion made 
in Ref. [ij- In the conclusions, we discuss the possibility 
of calculating the corrections semi-classically. 

Some words of caution. The separation between the 
NPP and the truncated series is a computational scheme. 
We do not claim that the NPP has a distinct physi- 
cal interpretation. In the past, there has been attempts 
|H S S Q to relate the difference between the average 
plaquette and its perturbative expansion to the gluon 
condensate used in calculations based on the Operator 
Product Expansion (OPE) in the continuum. The re- 
sults presented here can certainly be compared to the 
lattice calculations in these references, however our goal 
is to understand a particular lattice model studied with 
a particular perturbative scheme. We believe that an ac- 
curate parametrization of the NPP defined above as a 
function of /3, is a preliminary step that needs to be com- 
pleted before attempting to extract accurate numerical 
values of continuum quantities such as the gluon conden- 
sate. 



II. MINIMAL ERROR ESTIMATE 

In the following, we consider quantities Q which admit 
an asymptotic series of the form 

oo 

Q^^flfeA'^. (1) 

fc=0 

We assume that the leading growth of the coefficients can 
be parametrized as follows: 

au - CiC^r{k + Cs) (2) 

This type of behavior is justified generically in Ref. [l3| . 
In the case of lattice gauge theory, this type of behavior 
is seen explicitly in the limiting case of the one plaquette 
model discussed in Sec. |V] We define as the difference 
between the numerical value of Q and the series truncated 
at order k 

k 

Afc(A) = Q(A)-^azA' (3) 

1=0 

In addition, we will assume that for sufficiently small 
coupling, Afc is approximately given by the next order 
contribution 

Afc~ A'^'+^afc+i. (4) 

At fixed coupling A, |Afc| can then be minimized for 
k — k* with 



This expression has been obtained by using the lead- 
ing term of Sterling formula for the gamma function 
(V^z^-i/^e-^). In Eq. ©, (AlCal)"^ is the leading 
term and — C3 — (1/2) the first correction. It assumed 
that A is small enough to neglect the 0(A|C2|) terms 
which require the higher order terms of Sterling formula. 
Plugging Eq. (O in the expected error, gives the minimal 
error 

Miuk lAfel ~ V2^|Ci|(A|C2|)i/2-C3g-^ 
The exponential part of this formula is well-known [T^ . 

III. THE ANAHARMONIC OSCILLATOR 

A simple quantum mechanics example where the situ- 
ation described in section m is approximately realized is 
the anharmonic oscillator. The Hamiltonian reads 

H = + + Xx"^ . (7) 

We discuss the perturbative expansion of the ground 
state energy 

The leading asymptotic behavior of the coefficients has 
been calculated by Bender and Wu [T^ : 

afc~ (-l)'=+iv^3'=r(fc + l/2) (9) 

Using Eq. ||HJ), we obtain 

Miuk lAfcl ~ (VT2/7r)e-^ . (10) 

This is illustrated in Fig. ^ where the logarithm of the 
error is plotted versus 1/A. With this choice of variables, 
the minimal error is a represented by a straight line that 
is approximately tangent to all the curves representing 
the empirical error. 

A careful look at Fig. Qshows that the empirical curves 
go slightly below our estimated error. In other words, 
|Afe| is slightly smaller than Eq. This is illustrated 
in Fig. 121 One can see that the minimal error from Eq. 
(|10|l provide an envelope for the curves from Eq. Q at 
various order, but the empirical |Afe| are slightly smaller. 
One also sees that Eq. (0J become better as A becomes 
smaller. A more detailed study shows that the difference 
between the Eq. Q and |Afc| decreases as a power of A 
close to 1. The features of |Afe| seem related to the fact 
that the series has alternated signs. In the following, we 
will only consider same sign series and understanding the 
above feature is not crucial for the rest of the discussion. 

IV. THE DOUBLE- WELL 

A more intricate situation in encountered for the 
double-well potential. In shifted coordinates, the poten- 
tial reads 



k* ~ {X\C2\r' - C3 - (1/2) + 0{\\C2\) . (5) 



Viy) = {l/2)y'-gy^ + igy2)y* . (11) 
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FIG. 1: Natural logarithm of the absolute value of the differ- 
ence between the series and the numerical value for order 1 
to 15 for the anharmonic oscillator as a function of 1/A. As 
the order increases, the curves get darker. The thicker dark 
curve is In ((\/T2/7r)e~3>^) 
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FIG. 2: Natural logarithm of the absolute value of the differ- 
ence between the series and the numerical value for order 10 as 
in Fig. Q compared with the the estimate of Eq. @ at order 
10 (dashes). The thicker dark curve is In ((\/T2/7r)e~ 3a ) 
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FIG. 3: Natural logarithm of the absolute value of the differ- 
ence between the series and the numerical value for order 1 to 
15 (in g^) versus for the ground state of the double-well 
potential. As the order increases, the curves get darker. The 

thicker dark curve is ln( (g7r)"^/^e e^). The dash curve is 

In {^/W/Tvg~^e ^). 

If we expand the ground state in powers of g^, 

£;o^^afe5'\ (12) 

the leading asymptotic behavior of the coefficients reads 

El 

Ofc ~ -(3/7r)3'=r(fc + l) . (13) 

This implies 

Mink |Afc| - y/ej^g-^e"^ (14) 

Fig. O shows that this lower bound on the error is not 
reached. Rather, we see that the error is bounded by the 
one instanton contribution 

AEa = -(57r)-i/2e"^ . (15) 

This is a non-perturbative effect that at first sight is not 
related to the large order behavior of the perturbative 
series. However, some connection exists (see below). As 
illustrated in Fig. ^ a good approximation for the error 
at order k is 

Afe ~ ^{gn)-'^'c-^ + afc+ig^C^+i) . (16) 

Both terms are negative and there is no possible cancel- 
lation between them. 
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l/g^2 

FIG. 4: Natural logarithm of the absolute value of the dif- 
ference between the series and the numerical value for order 
10 (in g^) versus for the ground state of the double- well 
potential (solid line). The long dash curve is In (lang^^l). 

The dots represent In (| — ((/7r)^^''^e ^ + aiig^^\). 



Except for the overall normalization, Eq. H14|) is the 
square of the instanton effect. This is not a coincidence. 
In Ref. ^3 , it is explained that Eq. H13|l is an instanton- 
antiinstanton effect. It is possible to get rid of the one 
instanton effect by replacing the numerical value of the 
ground state by an average over the two lowest energy 
levels. This can be seen from the lowest order semi- 
classical formula [T6| : 

Eo ~ (l/2)-(g^)-i/2e-i^ , 
El ~ (l/2)-f 

One then recovers Eq. (|14|l as shown in Fig. [S| Note that 
the approximate doubling of the energy levels at small 
coupling is not seen in perturbation theory because one 
minimum of the potential goes to infinity when g goes to 
zero. 

As illustrated in Fig. a good approximation for the 
error at order k on the average of the two lowest energy 
states is 

Ar- ^ V6/^g-'e~^ + a^+ig^C^+i) . (17) 

The two terms are of opposite signs and there are possi- 
ble cancellations between them. This explains the spikes 
seen in Figs. Eland El 
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1/9-^2 

FIG. 5: Natural logarithm of the absolute value of the dif- 
ference between the series and the numerical value for order 
1 to 10 (in g'^) versus 1/g^ for the average between the two 
lowest energy states of the double-well potential. As the or- 
der increases, the curves get darker. The thicker dark curve 

is {gn)~^^^e ^ . The dash curve is \/&/iig~^e 
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l/g''2 

FIG. 6: Natural logarithm of the absolute value of the differ- 
ence between the series and the numerical value for order 10 
(in gf^) versus l/(?^ for the average of the two lowest states 
of the double- well potential (solid line). The dash curve is In 

(lvW^5"^e~^ -hang^^l) 
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V. A ONE PLAQUETTE GAUGE MODEL 



ONE PLAQUETTE 



Before discussing quenched QCD, we consider the sin- 
gle plaquette SU{2) gauge theory. After gauge fixing 
three of the four hnks, the partition function reads 



Z{(3) = / d[/c-^(i-5«'^i^-^) 



(18) 



Assuming /3 > 0, this can be rewritten as 

.2/3 



Z{p) = (2//3)3/2- / dtt^/\''^l-{t/2(3) (19) 



In Ref. ^3 1 it is shown that the integral can be expressed 
exactly as a converging expansion in /3^^: 



Z(/3) = (/37r)-3/22V2^A,(2/3)/3- 



fc=0 



with 



Ak{x) 



fc!(l/2-fc) ' 



(20) 



(21) 



Note that the peak of the integrand e~*i'^+^/^ becomes 
larger than 2(3 when k — 2(3 — 1/2, consequently the 
expansion is approximately truncated at that order. On 
the other hand, if we extend the range of integration from 
2(3 to +00, the integral in Eq. 1)21(1 becomes the gamma 
function and the coefficients grow factorially. We can 
determine the parameters Ci , C2 and C3 of the leading 
asymptotic behavior given in Eq. (0). We consider the 
sum without the prefactor (/?7r)~'^/^2^/^. From the study 
of ratios of successive coefficients we obtain C2 = 1/2 and 
C3 = 0. Using Stirling formula, we then obtain Ci — 1. 
From Eq. ((S)), the optimal order of truncation is 



k* ^2(3- 1/2 



(22) 



Note that this is exactly the same order as the order 
where the peak of the integrand moves outside of the 
range of integration in the exact Eq. H21() . From Eq. © 
and after restoring the prefactor, we obtain 



Mtuk lAfel ~ (21/2/71)/?- 



-2/3 



(23) 



Fig. [7| shows that this is a good approximation 

It is now possible to define precisely what we mean 
by "the non-perturbative part of Z(/9)". The perturba- 
tive part is the perturbative series (Eq. ((21|l with the 
argument of Ak replaced by 00) truncated at the integer 
closest to k* defined above, and denoted r(fc*) hereafter. 
The non-perturbative part consists of two parts. The 
first part, denoted R, consists in the remaining terms of 
orders r{k*) + 1 and higher in Eq. H2U|) . The second part, 
denoted T, is minus the integration tails for the first r(fc*) 
terms that we have added in order to get a perturbative 
series in 1/(3 with /3-independent coefficients. More ex- 
plicitly, 



ZW) 



^Pert. 



iP) + ZNPert.W) 



(24) 
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FIG. 7: Natural logarithm of the absolute value of the differ- 
ence between the series and the numerical value for order 1 
to 10 for the one plaquette integral as a function of (3. As the 
order increases, the curves get darker. The thicker dark curve 
is In ((2i/V^)e""'//3') 



with 



r(fc*) 



Zpert.W = (/37r)-3/22i/2 ^ A,(oo)/3-^ 



k=0 



and 



ZNPert.m = {Bi(3)-T{P)) 



(25) 



(26) 



i?(/3) = (/37r)-3/22i/2 ^ Afc(2/3)/3-'= (27) 

r(fc*) + l 



T{(3) = (/37r)-3/22i/2 
'^'^^ _,r(fc + l/2) 



(28) 



k=0 



fc!(l/2-fc) J,^ 



dte"V+'/^ 



The asymptotic behavior of T{(3) can be estimated us- 
ing section VI of Ref. p3| • We obtained that at leading 
order, the contribution of the tails up to order K is 



5Z{P,K) « Aife-2'3/3-i2^-3/2 



with 



r(^-i/2) 



(29) 



(30) 
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ONE PLAQUETTE 



ONE PLAQUETTE 





FIG. 8: In \R{f3) - T{l3)\ versus /3 (lowest, wavy line). The 
intermediate, almost straight curve is In {{2^^^ /Tv)P~^e~^'^). 
The dashed curve is \n{{2/-Kf^^l3-'^/^e-^'^) 



Writing (1 — e *) as e *(e* — 1) and expanding e* we 
obtain 

g r([_i/2) ^ dtt-'/2^^ _ g-t) ^ 2^1/2 . (31) 
1=1 ■ 

For large K, 

E^^^^^- rrftt-/^(l-e-*), (32) 
1=1 -^0 



we obtain the leading behavior 



A 



K 



(33) 



K 



Replacing K by its optimal value in Eq. (|22|1 . we con- 
clude that for large (3 



Tip) ~ (2/^)3/2/?-3/2e-2/3 



(34) 



It is interesting to notice that T{f3) is larger by a factor 
/3^/2 than the estimated error given in Eq. This 
means that we should also have 



i?(/3) ~ (2/^)3/2/3-3/2e-2/3 



(35) 



so that the two leading order contributions cancel. The 
correctness of this argument is illustrated in Figs. Inland 

m 

In conclusion, we have learned that in this simple lat- 
tice model with a compact gauge group, the non per- 
turbative part of the integral has two pieces. One piece 



FIG. 9: In |i?(/3) - r(/3) | + 2(3 versus ln(/3) (wavy curve in the 
lower part of the graph). The thicker line is In {{2^^^ /Tv)/f3^). 
The dashed line is \ii{{2/TTf^^ (3-^^^). In \R{3)\+2(3 (points) 
and In |r(/3)| -I- 2/3 (line) are inter weaved between the two 
straight lines. 



comes from the higher orders, is positive and dominates 
at sufficiently large (3. The other comes from the added 
tails of integration, is negative and dominates at suffi- 
ciently small (3. It is possible that this pattern persists 
for the usual lattice gauge models on cubic lattices. 



VI. QUENCHED QCD 

We now proceed to introduce lattice quenched QCD, 
the main model discussed in this article. With standard 
notations, the partition function is 



Z = 



n / dUie 



(36) 



with P — 2N I . For symmetric finite lattices with 
sites and periodic boundary conditions, the number of 
plaquettes is 



TVp = L°D{D - l)/2 
Using the free energy density 

/ = -(l/AAp)lnZ, 
we define the average plaquette 

P = df/dp , 

and its perturbative expansion 



(37) 
(38) 
(39) 
(40) 
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QUENCHED QCD 
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FIG. 10: Natural logarithm of the absolute value of the dif- 
ference between the series and the numerical value for order 
1 to 10 for quenched QCD as a function of /3. As the order 
increases, the curves get darker. The approximately straight 
dark curve is In (6.5 x 10« x (^)204/i2ig-(i6.V33),3) 



QUENCHED QCD 




FIG. 11: Natural logarithm of the absolute value of the dif- 
ference between the series and the numerical value for or- 
der 5 for quenched QCD as a function of /3 (dots). The 
dash line is hi{\ae/ f3^\). The solid curve is In (6.5 x 10** x 

(^)204/121g-(167rV33)/3^ 



In the following, we used the first three analytical values 
from Ref. 18], and the coefficients from 4 to 10 from Ref. 
0. The coefficients of Ref. up to order 16 are given 
in a figure and will be used to check extrapolations. At 
large /3, high precision data is necessary and we have used 
some values of the average plaquette from Refs. |0, . 

A simple guess is that the envelope of the accuracy 
curves is given by a power of some renormalization group 
invariant scale. For instance, one could test the idea that 
it is proportional to the fourth power of the two-loop 
renormalization group invariant scale 

Mink lAfcl ^ C(/3)204/i2ie-(i6.V33),3 (4^) 

Fig. II 1)1 shows that up to order 10, this provides a reason- 
able envelope in the region 5.5 < /? < 6. The constant 
C has been fixed to 6.5 x 10^ by finding a plateau in 
|Aio|(/3)-^°''/^2^e(i6'^'/33)'3. As the curves leave the con- 
jectured envelope, they become more fiat. For instance 
at order 8 and for 6 < /? < 7, a reasonable fit can be 
obtained 6] with the square of the perturbative renor- 
malization invariant scale. However, this does not seem 
to characterize the asymptotic behavior. Note also that, 
at values of /3 close to 6, the empirical error is signifi- 
cantly larger than the next order contribution. This is 
illustrated in Fig. ^]at order 5. 

Another simple guess is that the envelope is propor- 
tional to a power of the lattice spacing a expressed in 
units of tq = 0.5 fm. For the interval 5.7 < (3 < 6.92, the 
following empirical power series [To( was obtained from 



the force method. 

ln(a/ro) = -1.6804 - 1.7331 (/3 - 6) 

+ 0.7849 (/3 - 6)2 - 0.4428 (/3 - 6)^ . 

It has been suggested that the envelope is proportional 
to 0i S J however this is hard to establish without 
the knowledge of the higher orders. On the other hand, 
the accuracy curve at order 10 can be fitted very well 
with a? but again it does not seem to characterize the 
asymptotic behavior since the successive orders are still 
quite separated. This is illustrated in Fig. [T^ 

The question is now, if the higher orders can be calcu- 
lated will the envelope stay approximately straight as in 
Fig. ^]or curl up as in Fig. II 21 ' These two alternatives 
need to be compared systematically. As we will see in the 
next section, this requires the knowledge of coefficients 
of higher orders than those actually calculated. 

VII. LARGE ORDER EXTRAPOLATIONS 

In this section we consider two large order extrapo- 
lations that can be built out of the first 10 terms used 
above. The first extrapolation is based on the assump- 
tion justified in Ref. 0] that dPjdfi has a logarithmic 
singularity in the complex /3 plane. Integrating, we ob- 
tain 

^ ^ C(Li2(/3-V(/3-' + ^F)) + h.c , (42) 

fc=i 



QUENCHED QCD 
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FIG. 12: Natural logarithm of the absolute value of the dif- 
ference between the series and the numerical value for order 
1 to 10 for quenched QCD as a function of /3. As the or- 
der increases, the curves get darker. The short dash curve is 
ln(0.14 (a/rof), The long dash curve is ln(1.5 (a/ro)*). 

with 

Li2(a;) =^xVfc^ (43) 
fc=i 

In Ref. [3 we argued that 0.001 < T < 0.01 and for this 
range of values, the low order coefficients depend very 
little on r. For this reason, it is difficult to fit its value 
given the accuracy of the coefficients. We have taken 
the intermediate value T = 0.003 which is compatible 
with everything we know and determined C = 0.0654 
and f3rn=^-7&7 using the known values of ag and oio on 
a 24"^ lattice 2]. The estimated errors on these coeffi- 
cients and their impact on the determination of C and 
/3m are discussed at length in the next section. Other 
choices of F with this order of magnitude do not affect 
our conclusions. The numerical values are given in Table 
I. Except for the first term, the agreement is very good. 
The fact that such a good agreement can be reached 
by tuning two parameters begs for a diagrammatic ex- 
planation! We also obtained oig — 7.7 10^. This is 
quite close to what we found from Fig. 1 in Ref. 0: 
ai6 ~ 0.00027 X ~ 7.6 10^ Note however, that this 
value of oig was obtained on a 8'* lattice and the finite 
volume effects should be significant. The perturbative 
series obtained from Eq. H42|l has finite radius of con- 
vergence 3, 8J. This seems to contradict the common 
wisdom [21| that if we ignore non-perturbative effects 
(responsible for the fact that P takes different limits 
when — > 0^), we should run into much more serious 
problems. Consequently, we believe that Eq. H42|) can 
only be a good model for the low orders. 
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m 


bm. 


bm 


1 


2 


0.7567 


2 


1.2208 


1.094 


3 


2.9621 


2.811 


4 


9.417 


9.138 


5 


34.39 


33.79 


6 


136.8 


135.5 


7 


577.4 


575.1 


8 


2545 


2541 


9 


11590 


11590 


10 


54160 


54160 



TABLE I: bm'- regular coefficients; bm'- extrapolated coeffi- 
cients for r = 0.003. 

Models based on infra-red renormalons 0, 0, 0| have 
a better chance to reproduce the large order behavior of 
the series. We assume 

Pf=iK r dte-/^' (1 - t 33/167r2)-i-^ (44) 
Jti 

with 

/5 = /3(l + di//3+...) (45) 

In Refs. Ulil, the value x = 204/121 was used. Note 
that ti = corresponds to the UV cutoff (we use that 
value later) and ^2 — 167r^/33 corresponds to the Landau 
pole. Expanding (1 — t 33/167r^)~^~^ in powers of t and 
extending the integration range to c», we find that at 
leading order 

PW - E ' (46) 

m—l 

with 

bk = Ci (33/167r2)'= T[k + x] . (47) 

Setting ^2 and higher order coefficients to zero, x — 
204/121 and using the known values of ag and aio , 
we obtain Ci = 0.219 and di = —3.82. The low order 
predictions are less good than for the previous method, 
for instance, ag = 120.4. At larger order we obtained 
ai6 = 7.4 10® which is close to 7.6 10®. Note that 
g47r"(3.82)/33 ^ gj jg significantly larger than the usual pl | 
factor 28.81 found for the conversion to the MS scheme. 

The two extrapolations models are compared in Fig. 
1131 The two models yield similar coefficients up to order 
20. After that, the integral model has the logarithm of 
its coefficients growing faster than linear. 

Using Eq. ©, for the expansion, we obtain 

Mzuk \Ak\ ~ 3.5(^)204/i2i-i/2e-(i6.V33)/3 (48) 

Shifting to /3 using Eq. and neglecting (3~^ correc- 
tions 

Mlflk lAfel ~ 3.1 108(/5)2O4/121-l/2g-(16.V33),3 (49) 
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FIG. 13: ln( b^) for the dilogarithm model (solid line) and 
the integral model (dashes). The dots up to order 10 are the 
known values. 



Except for the -1/2 in the exponent, this is proportional 
to the two- loop renormalization group invariant. The 
situation is reminiscent of what we encountered in Sec. 
IVl where the error was smaller by a factor /?^^/^ than 
naively expected. Another possibility would be to take 
X = 204/121 -|- 1/2 if we insist on having a quantity that 
is scheme independent. 

The study of the asymptotic behavior of the empirical 
f]'-^ expansion between orders 40 and 60 indicates C2 — 
0.203 quite close to 33/(167r2) ~ 0.209. The other values 
obtained are quite different from those of the coefficients 
of the expansion: Ci ~ 990 and C3 ~ 4.4. With 
these values, the optimal order for /3 = 6 is 25. 

The accuracy curves for the two models are compared 
in Figs. El and El We have used the known values 
for the first 10 coefficients. A careful examination shows 
that the figures hardly differ up to order 25. From these 
two figures, it is plausible that the envelope of the true 
series is close to 0.65 (a/ro)^. This value of 0.65 does not 
correspond to a best fit procedure. If we had forced the 
{a/ro)^ behavior and fitted the proportionality constant 
at order 26, we would have obtained 0.79 for the diloga- 
rithm model and 0.71 for the integral model. If we had 
plotted the corresponding curves in Figsll4l and ll5l they 
would hardly be visible on these black and white figures. 
By taking a slightly smaller value 0.65, we made the dash 
curve more easy to see. We could also have plotted the 
curve corresponding 0.4 (a/rp)* that can be inferred from 
the graphs and the estimate of the gluon condensate of 
Ref. 01 . Such a curve would be -0.49 below the one we 
drew. 

The tadpole improvement 24J could be another method 



FIG. 14: Natural logarithm of the absolute value of the dif- 
ference between the series and the numerical value for order 
1 to 30 for quenched QCD with the dilogarithm model. As 
the order increases, the curves get darker. The long dash 
curve is ln(0.65 (a/ro)*). The solid curve is In (3.1 x 10* x 
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FIG. 15: Natural logarithm of the absolute value of the dif- 
ference between the series and the numerical value for or- 
der 1 to 30 for quenched QCD with the integral model. As 
the order increases, the curves get darker. The long dash 
curve is ln(0.65 (a/ro)*). The solid curve is In (3.1 x 10* x 



204/121-1/2 -(167r^/33)/3\ 
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to study the asymptotic envelope of the accuracy curves. 
One defines the new series 



P 



K 

m=0 



K 



-K-l^ 



m=0 



with a new expansion parameter 



1 



(50) 



(51) 



The consistent sets of numerical values 0, 0, ll^ , show 
that the beginning of the series converges much faster. 
This is clear up to order 8. However, for larger orders, the 
dependence on the higher order coefficients of the origi- 
nal series due to powers of 1/(1 — J2m=o ^m/'^™) lo'^ 
orders in /3jj makes the tadpole improved series almost as 
slow as the original series. In extracting the coefficient 
of a^, P. Rakow used a more sophisticated method where 
1/(1 — J2m=o bmfi~^) is replaced by the improved series 
with a value of (3r determined numerically |2^. This 
procedure gives better perturbative estimates which ex- 
plains that the constant of proportionality 0.4 is smaller 
than 0.65 as discussed above. 

Instead of assuming the (a/ro)'* behavior, it is also pos- 
sible to extract the power from linear fits of ln(| A^j) ver- 
sus ln(a/ro). For fc < 26 the two extrapolations give sim- 
ilar fits near (3 — &. For instance, A26, is 0.56(a/ro)^' ''^ 
in the dilogarithm model and 0.63(a/ro)'^'^*' for the inte- 
gral model. For larger k, the constant of proportionality 
and the exponent keep increasing slowly for the diloga- 
rithm model and the procedure becomes meaningless for 
the integral model. We now discuss into more detail the 
stability of this fitting procedure. 



VIII. UNCERTAINTIES ON 
EXTRAPOLATIONS 

We now discuss the effects of the uncertainties of the 
known coefficients on the extrapolations introduced in 
the previous section. The determination of the unknown 
parameters in both models depends only on ag and aiQ. 
The errors on these coefficients are due to the finite vol- 
ume and the numerical implementation of the stochastic 
perturbation theory. According to Ref. 0, from which 
we took the values of og and oio for a 24** lattice, the 
statistical errors on ag and aio are close to one percent 
while the errors due to the finite volume are less than 
half a percent for the calculations on a 24^ lattice. 

We have studied the effects of small changes of the 
following form 



ai —f aj(l + 5i), 



(52) 



for 1 = 9 and 10 on the determination of the parameters 
A (amplitude) and S (slope in a log-log plot) appearing 
in the approximate parametrization 



which was our best estimate of the NPP of the plaquette 
in the previous section. With the sign convention for A^. 
in Eq. 0,1. e., the numerical value minus the series 
at order fc, A^ > for k < 10. If the coefficients of the 
series are all positive and end up growing factorially, it is 
clear that for a given /3, there is some order for which the 
truncated series becomes larger than the numerical value 
and Afc becomes negative. This explains the downward 
spikes that are clearly visible for orders 29 and 30 of the 
integral model on Fig. El On the other hand, for the 
dilogarithm model, A^ remains positive as k increases, 
over the whole P interval for which A^ is larger than the 
errors on the numerical values of P. 

When the value of the series becomes close to the nu- 
merical value, the errors on the Monte Carlo estimates 
of P become important. The numerical values that we 
have used for the fits of A and S have estimated errors 
01 smaller than 5 x 10~^. Consequently, the details of 
Figs. 1 1 41 and [T51 b elow -10 on the y axis are meaningless. 

We have performed fits of A and S using ln(|A26|) 
versus ln(a/ro) for 18 values of (3 between 5.68 and 6.10. 
The quality of the fit can be assessed from 



1/2 



Xf 



(l/18)£(A/ln(|A26(A)l))' 



(54) 



1=1 



with 



A = ln(|A26(A)|) - In(^) - 51n(a(ft)/ro) . (55) 

Xf give an estimate of the typical relative errors obtained 
with the fit. 

It is also useful to assess how well a particular choice 
of parameters reproduces the known values of the known 
coefficients of order 8 and less. For this purpose we have 
defined 



1/2 



Xa = I (l/4)^[(aj -aj[ag,aio])/ajf 

j=5 



(56) 



A26 ^ A{a/roY 



(53) 



with [09,010] understood as the value of aj obtained 
from the models after the unknown coefficients have been 
fixed using ag and oig. We started at order 5 because the 
differences on lower orders are quite large for the integral 
model. 

It is interesting to notice that in both models, the ex- 
trapolated coefficients depend linearly on a overall scale 
factor called C in the dilogarithm model and Ci in the 
integral model. These scale factors cancel if we consider 
the ratio of the two equations determining the unknown 
parameters. In other words, the parameters (3m and di 
depend only on ag/oig. Consequently if ^g = 5io — 6, 
ag /oio is unchanged and the only effect of the change is 
to rescale all the extrapolated coefficients by a common 
factor 1 -I- ^. This type of modification generates a lin- 
ear response which is the same at all order and has quite 
controllable effects. On the other hand, if ^g 7^ Siq, the 
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59 


Sio 


c 


/3m 


A 


5 


Xf 


Xa 








0.0654 


5.79 


0.56 


3.72 


0.004 


0.006 


0.01 


0.01 


0.0661 


5.79 


0.60 


3.81 


0.005 


0.006 


-0.01 


-0.01 


0.0647 


5.79 


0.52 


3.65 


0.003 


0.015 


0.01 


-0.01 


0.0791 


5.67 


0.40 


3.31 


0.002 


0.048 


-0.01 


0.01 


0.0541 


5.90 


1.46 


4.72 


0.017 


0.056 


-0.005 


0.005 


0.0595 


5.85 


0.77 


4.07 


0.007 


0.031 



TABLE 11: Values of C, f3m, A, S, Xfy a-nd Xa defined in 
the text for the dilogarithm model, corresponding to various 
choices of Sg and Sio. All the fits refer to a comparison with 
the 26th order. 



59 


5io 


Ci 


~di 


A 


5 


Xf 


Xa 








0.219 


3.82 


0.63 


3.90 


0.004 


0.08 


0.01 


0.01 


0.221 


3.82 


0.70 


4.01 


0.006 


0.07 


-0.01 


-0.01 


0.216 


3.82 


0.57 


3.81 


0.003 


0.08 


0.01 


-0.01 


0.263 


3.72 


0.42 


3.44 


0.002 


0.03 


-0.01 


0.01 


0.182 


3.92 


2.29 


5.24 


0.024 


0.12 


-0.005 


0.005 


0.199 


3.87 


0.94 


4.33 


0.009 


0.10 



TABLE III: Values of Ci, -di, A, S, Xf, and Xa defined 
in the text for the integral model, corresponding to various 
choices of 5^ and Sio- All the fits refer to a comparison with 
the 26th order. 



parameters (3rn and di are modified and this has more 
drastic effects on the extrapolated coefficients. 

In order to give a quantitative idea, we have con- 
sidered the four cases corresponding to the condition 
ISg] = \Sio\ = 0.01. The results are shown in Tables 
II and III. The choice = —0.01 and 5io = 0.01 has a 
Xf significantly larger than the other choices. For this 
choice, the two models start differing significantly at or- 
der lower than 26 and then develop negative spikes at 
larger order. It might have been a better choice to quote 
fits at lower order where there is a better agreement be- 
tween the two models and smaller Xf- For instance, at 
k — 22, we obtain A — 0.69 and S — 3.90 with Xa — 
0.007 for the dilogarithm model and A — 0.61 and S = 
3.77 with Xa ~ 0.005 for the integral model. As the fit 
for fc = 26 did not look reasonable for this choice of Si, 
we also added a line with the same signs but a lower 
magnitude. 

In summary, we obtained values of S between 3.3 and 
4.3 for fits with a reasonable Xf- The boundary of this 
interval corresponds to fits that have a larger Xa- These 
results seem to favor the idea that the NPP scales like a'* 
over the idea that the NPP scales like a^. However, no 
strong conclusions can be drawn because, the extrapola- 
tion models are empirical. The only indication of their 
approximate validity is that they provide consistent re- 
sults up to some reasonably large order. It should also be 
pointed out that there is no universal agreement 0, |23| 
about what should be done about the Landau pole in Eq. 
(|44|l . In order to draw strong conclusions, one would need 
either, numerical values of the coefficients up to signifi- 



cantly larger order than available now, or extrapolation 
models supported by a detailed analysis of graphs having 
significant contributions in lattice perturbation theory. 



IX. PARAMETRIZATION OF THE FORCE 
SCALE 

In the previous section, we considered the possibility 
of having the minimal error proportional to the fourth 
power of the two loop renormalization group invariant 
scale. What is attractive about this possibility is that it 
has the generic form A e~'-^/'^. However, the extrapo- 
lations suggest that it is not a very accurate parametriza- 
tion of the envelope of the accuracy curves and that a 
(a/ro)'* parametrization seems more promising. In this 
section, we consider simple parametrizations of a. 

First we study the derivative of ln(a) with respect to (i 
with the first two universal terms subtracted. It should 
be noted that the derivative emphasizes the cubic terms 
in the expansion of ln(a) and the difference between the 
two expansions from Refs. and [Toj | is more pro- 
nounced than without the derivative (see Fig. I16|l . We 
propose the following parametrization. 

d\n{a/rn)/d(3 = -(47rV33) + (51/121)/3-i + AiQ-'^^'^ . 

(57) 

This form was originally motivated by our idea of obtain- 
ing a NPP that could be calculable semi-classically and 
was found later in good agreement with an independent 
argument ^ (see below). We have also tried power 
parametrizations of the form Bi(3^^ and found that it 
requires a large value of B2 (close to 16) not very sta- 
ble under small changes of the fiting interval. Using the 
Taylor expansion of Ref. |lO| to produce a set of numer- 
ical values between 5.9 and 6.3, subtracting the first two 
terms of the r. h. s. and taking the log, we can determine 
numerically the parameters Ai and A2 from a linear fit. 
The result is Ai = -1.35 10^ and A2 = 2.82. For these 
values the exponential corrections become larger than the 
two universal terms for (3 < 5.8. Integrating and fitting 
the constant of integration in the same interval, we found 



ln(a/ro) = A3-(47rV33)/3+(51/121)ln(/3)-(AiM2)e-^^^ 

(58) 

with A3 — 4.5276. This constant was obtained by plot- 
ting the 1. h. s. minus all the other terms of the r. h. s. 
versus p. A nice plateau appears between /? = 5.9 and 
6.3. The extremal values in this interval are 4.5272 and 
4.5282. 

Note that the value of A2 seems consistent with the 
idea of using Op^^j corrections for this quantity. The 
symbol Qpert refers to the one-loop or two-loop expres- 
sion which in the short f3 interval considered here can 
hardly be distinguished from each other. For this reason, 
we have not introduced any power correction in the last 
term of Eq. (|57|) . The assumption of Cpg^j corrections 
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dln(a/r_0) /d/3 
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FIG. 16: d\vL{a/ro)/d(i using Ref. O (thick dashes), ||3 
(smaU dashes) and Eq. 1571 (solid hne). 

fixes A2 — 87r^/33 ~ 2.4 which is close to the value 2.82 
obtained above. 

It is possible to use Eq. (|55|) to predict ln(a/r'o) at 
large (3. For instance, at /? = 7.5, we obtain -3.59 (for 
-3.63 in Ref. I^) and -4.74 at /3 = 8.5 (for -4.81 in 
Ref. [23 ). It is also possible to obtain the lattice scale 
Ai from the constant of integration A3, namely A^, = 
exp(-yl3)/ro ~ 4.4 MeV. 

X. CONCLUSIONS 

We have considered the definition of the non- 
perturbative part of a quantity as the difference between 
its numerical value and the perturbative series truncated 
by dropping the order of minimal contribution (which 
is coupling dependent) and the higher orders. For the 
anharmonic oscillator, the double-well potential and the 
single plaquette gauge theory, the non-perturbative part 
can be parametrized as A e~^/^ and the constants 
A, B and C can be calculated analytically. 

For lattice QCD in the quenched approximation, the 
perturbative series for the average plaquette is dominated 
at low order by a complex singularity and at this point, 
the asymptotic behavior can only be reached by using ex- 



trapolations. We have considered two extrapolations that 
provide a consistent description of the series up to order 
20-25. More work is needed to understand these extrap- 
olations better. A diagrammatic interpretation of the 
dilogarithm formula should be possible. The similarities 
between the one plaquette integral and the renormalon 
motivated integral suggest that it might be possible to 
derive the latter from an effective action for a single pla- 
quette, obtained from the large volume partition function 
after integrating over all the other links. Ultimately, the 
two approximations should be put together using disper- 
sion relations in the complex (3 plane. 

The two extrapolations favor the idea that the non- 
perturbative part is a power of the lattice spacing cal- 
culated using the force method and that the power is 
close to 4. We found a parametrization of the force scale 
as the two loop universal terms with exponential correc- 
tions consistent with flpg^^ corrections. These corrections 
become quite important when (3 is near or below 6. 

For large /9, we can treat the exponential corrections 
as small quantities and expand in power of these small 
quantities. The NPP of the plaquette is then written as 
superposition of terms of the form Aj (3^^ e~'-^^^. If the 
constants Bj and Cj can be determined from general ar- 
guments, we only need to determine the amplitudes Aj, 
starting from those with smallest Cj. Calculating more 
terms of the perturbative series and calculating more ac- 
curate numerical values of the average plaquette at large 
(3, would help in the empirical determination of these con- 
stants. As /3 becomes large, the lattice spacing becomes 
small and the NPP is small. It is plausible that the am- 
plitudes could be calculated by semi-classical methods in 
the continuum. This is the most challenging part of the 
program. A method that comes to mind is the instanton 
calculus. However, the infrared sensitivity of the pro- 
cedure is notorious and phenomenological input such as 
the gluon condensate is often invoked to produce realistic 
formulas p9|. Having accurate empirical formulas for a 
well defined model (quenched) lattice QCD would help 
to revisit this question from scratch. 
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